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1  2 

In  recant  works  by  the  authors  '  the  problem  of  charge-carrier 

transport  in  a  semiconductor  in  the  presence  of  an  applied  electric 

field  was  examined.  A  new  nonlinear  differential  equation  for  the 

distribution  function  of  charge  carriers  was  derived  incorporating 

acoustic-strain  interactions  and  solved  exactly.  In  the  present  work 

this  analysis  is  extended  to  include  equivalent  intervalley  scattering 

3  4 

which  is  significant  to  transport  processes  in  Silicon.  '  The 
analysis  considers  only  single  phonon  interactions.  Furthermore, 
charge-carrier  energy  is  sufficiently  small  so  that  non-equivalent 
intervalley  transitions  are  omitted.  Again  a  differential  equation  is 
derived  for  the  distribution  of  charge  carriers  as  function  of  energy. 
However,  in  the  present  case,  this  equation  includes  ''difference”  terms 
stemming  from  intervalley  scattering.  An  approximate  solution  to  this 
e^idation  is  obtained  in  the  following  manner.  As  a  first  step,  a  group 
of  difference  terms  is  approximated  by  a  differential  form  whose 
components  are  motivated  through  matching  with  correct  asymptotic 
solutions.  The  resulting  nonlinear  equation  is  solved  exactly  thus 
yielding  an  approximate  solution  to  the  original  equation. 

In  the  second  step  this  approximate  solution  is  employed  in 
constructing  a  solution  to  the  full  nonlinear  difference-differential 
equation.  The  resulting  distribution  appears  as  a  piecewise  continuous 
function  separated  at  the  intervalley  phonon  energy.  The  solution  so 
obtained  shows  very  good  agreement  with  results  generated  by  a  Monte 
Carlo  analysis  of  the  problem.  The  new  analytic  distribution  is  used 
to  obtain  a  closed  expressions  for  drift  velocity  as  function  of 
applied  electric  field,  charge-carrier  density  and  temperature. 


Comparison  of  these  results  with  experimental  values  is  found  to  give 
very  good  agreement  for  n-type  Silicon  at  room  temperature  and  fields 
up  to  10^  V/cm. 
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II.  ANALYSIS 
A.  Starting  Equations 

In  this  work  we  consider  the  following  generalized  Uehling 

Uhlenbeck^'®  equation  for  the  distribution  function,  f(k,t), 

3f  .  eE,  df  '  ^  . 

3t  *  3k  ”  ^SA^^^ 

A 

Where  Jg^  represents  charge  carrier  -  strain  acoustic  interactions  and 
represents  equivalent  intervalley  interactions.  Electric  field  and 
particle  momentum  are  denoted  by  E  and  to  k,  respectively.  Both 

A# 

collision  terms  are  of  the  form 

2  [f'(l-af)  -  f(l-af)  (2) 

Where  f  s  f(k')/  etc.  and  represents  the  scattering  from  the 

state  k'  to  k  due  to  the  mechanism  a.  The  bookkeeping  parameter  a  has 

the  value  1  and  may  be  set  equal  to  zero  in  the  classical  domain.  In 

writing  (2)  we  have  assumed  the  normalization 

2f(k)  -  N  (3a) 

k  '' 

In  the  continuum  limit  this  relation  becomes 

— ^  /  f(k)dk  -  N  (3b) 

where  V  is  crystal  volume. 

In  this  analysis  we  will  concentrate  on  the  transport  of  charge 
carriers  in  Silicon.  The  possible  phonon-charge  carrier  interactions 
for  this  case  are  shown  in  Fig.  1.  It  is  assumed  that  charge-carrier 
energy  is  sufficiently  small  {i  1  eV)  thereby  permitting  us  to  neglect 
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non-equivalent  intervalley  scattering.^  In  this  diagram  "S”  denotes 
"strain”  (often  refered  to  as  "deformation  potential") ,  0  and  A 
represent  "optical"  and  "acoustic"  phonons,  respectively.  In  Silicon 
there  are  two  possible  intervalley  (IV)  transitions  denoted  by  g 
(parallel  valleys)  and  f  (perpendicular  valleys) .  The  terms  T  and  L 
refer  to  transverse  and  longitudinal  phonons,  respectively. 

To  determine  the  relative  importance  of  scattering  mechanisms  we 
define  the  total  scattering  rate 


x(a)  _  y  c(a) 

■  r,  k' 


(4) 


These  rates  are  calculated  in  Appendix  A  and  are  sketched  in  Fig.  2 

relevant  to  Silicon  at  300°K.  It  is  evident  from  these  figures  that 

the  interactions  g  (L0±)  and  (SA±)  dominate  thereby  justifying  the  RHS 

of  (1).  In  this  notation  (+)  refers  to  phonon  absorption  and  (-)  to 

phonon  emission.  It  should  be  noted  that  the  SO  interaction  listed  in 

Fig.  1  does  not  appear  in  the  rate  calculations  of  Fig.  2.  This  is 

attributed  to  symmetry  properties  of  the  crystal  which  in  turn  disallow 

4  8 

this  transition.  ' 

B.  Reduction  of  Kinetic  Equation 

To  account  for  anisotropy  of  the  distribution  function  we 

9 

introduce  the  expansion 


f(k,t)  -  2  ft(«,t)  P^(u) 


(5a) 


where 


,  -  liV 

®  “  2m 


H  2  E*k 


(5b) 


6 


J 


Hatted  variables  are  unit  vectors  and  m  is  effective  mass.  Note  that 
in  writing  (5b)  the  constant  energy  surface  has  been  approximated  by  a 
sphere.  Going  to  equilibrium  and  keeping  only  the  first  two  terms  in 
(3)  we  obtain 

H  a  It  («ti)  +  («) 

^  V  ItlG 


.  2 

For  we  obtain 
SA 


-  ^Sa(^o)  -  ‘1 


(7) 


where  the  relaxation  time  Tg^  is  as  implied  and 
g(f  )  s  k^T  ^  +  f  (1  -  af  ) 


(7a) 


Acoustic  mean  free  path,  t,  is  given  by 

t  =  _ 

m^Ej  kgT 


1,3 


(7b) 


In  this  expression  p  is  crystal  mass  density,  u  is  acoustic  phonon 
velocity,  T  is  crystal  temperature  and  Ej^  is  the  deformation  potential 
constant.  In  obtaining  (7)  it  was  assumed  that 

'fiuq  <<  kgT  (8) 


where  q  is  phonon  wavevector.  For  Jj^(f)  we  obtain  (see  Appendix  B) 


-  (t„  .  4ti)(l  -  af,)(njv  *  I  +  fl! 


(9) 
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i 


where  the  sum  is  over  upper  and  lower  terms  in  the  stimmand  and  are 
associated  with  phonon  emission  and  absorption.  Furthermore 


^  rti 


(9a) 


Here  we  have  introduced 

9(x)  »  1,  X  >  0 

0(x)  -  0,  X  <  0 


(9b) 


The  displaced  distribution  is  given  by 
f^(«)  »  ± 


(9c) 


and  the  length  is  written  for 


2pu 


'IV 


D^m^ 


(9d) 


where  D  is  the  intervalley  interaction  constant  (in  units  of 
3  4 

energy/ length) .  '  Intervalley  phonon  frequency  is  written  The 


intervalley  phonon  density  is  taken  to  be 

1 


n 


IV 


exp 


IV 


KbT 


(9e) 


-  1 


Note  that  (9)  has  the  structure 


Jiv(f)  =  Jiv(f^)  -  ^  f. 


(10) 


3 


I 


.  I  aj  HtoOiv  *  T  *  5) 


w  2  2  *  ^  ^ 


'o<"iv  *  ?  ♦  i  *  »  foU 


(10a) 


where  again  is  as  implied.  Inserting  (6,7,10)  in  (1)  and  equating 
Legendre  coefficients  we  obtain 


i  (if  \  m  ^  3  r  X  i  ^  \ 

3  (3<  '*^1^  eEi  k_T  d£  9(fo) 

O 


(lla) 


*  h  -  'o'"  J  *  5  *  "  1 


(11b) 


at 


-=—  f  -  I  — 
eEt  ^1  “ 


eE  “±^"lV  1 


-r  1 


t  *  ^  ^o)  ^ 


Eliminating  f^  from  (lla,b)  and  introducing  the  nondimensional  energy 
X  5  «/kgT  (12) 


gives  the  desired  nonlinear  difference-differential  equation 

±  ■^2*2' 


(♦fo')'  [x^  5(f^)]'  +  2  i  i  ?) 


(13) 


^0^"lV  2  2  *  *  ^O^  ^ 


where  prime  denotes  differentiated  with  respect  to  x. 

♦  H 


Furthermore 


sx' 


fi  S  r(nj^  I  ^  I  *  a  ^x(x±A)  9(x±A) 
s  rt  >/x(x±A)  9(x±A) 


(14) 
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In  the  proceeding  we  have  introduced  the  dimensionless  constants 

"ftw. 

r  3  i/6„.  .  A  5 


’IV 


IV 


,  a  ^ 

2mu' 


and  dimensionless  electric  field 


3,8 


-  (eEM^  =  ,E_  2 


6mu  kgT 


E 


(15b) 


Note  in  particular  that  setting  r  *  0  in  (13)  has  the  effect  of 
shutting  off  intervalley  scattering. 

With  these  parameters  given  by  (15),  solving  (11b)  for  f^  and 
inserting  the  result  into  (5a)  gives 

f(X,4)  »  [1  -  4  ^  f^(x)  (16) 

The  noramlization  (3b)  assumes  the  following  form  when  written  in  terms 

of  the  nondimensional  energy  x: 

/f  (x)Vxdx  -  Vir  A  (17) 

2 

where  A  is  the  quantum  degenerary  parameter^*^ 

A  =  n^\^^  (17a) 


and  is  the  thermal  deBroglie  wavelength 


mkgT 


(17b) 


and  n  denotes  charge  carrier  denisties.  Some  typical  values  of  these 

4 

parameters  are  displayed  in  Table  1. 

C.  Limiting  Properties 

We  wish  to  consider  solution  to  (13)  in  three  limiting  cases: 

(a)  r  =  0 


10 


Thus  has  the  effect  of  turning  off  intervalley  interactions  (i.e.,  L 
-*<*>).  In  this  event  (13)  reduces  to 


11 


I 


■ ; 

Values  of  Physical  Constants 
Used  in  This  Analysis 


material 

n-type  Si 

T 

300°K 

A 

2.40 

"iv 

0.100 

7 

85.5 

r 

1.65 

T 


497  V/CM 


(18) 


w 


(sxf^)  •  +  ]  •  -  0 


whose  solution  was  studied  in  detail  in  reference  1  and  is  given  by 


a  ^ 


(19) 


where  B  is  a  normalization  constant. 

■Cfe)_s_= 

This  limit  corresponds  to  turning  off  the  applied  electric  field. 
There  results 

g(f^)  5  2  C^(fQ)  »  0  (20) 


where 


Si<tol  =  Ci  t£,(niv  *  5  *  5) 


(20a) 


A  physically  relevant  solution  to  (20)  is  given  by  the  Fermi-Dirac 
distr  bution 

1 


a+Be 


(21) 


Note  in  particular  that 


(21a) 


(C)  A-»0 


This  limit  corresponds  to  <<  kgT. 

"iv  -  i  -  5  * 


In  this  limit  we  find 

(22a) 


^o  “  ^o  *  ^^o’ 


(22b) 


13 


in  which  case 


^  i  [g(fo)]'  +  O(A^) 


Furthermore,  with  8(x  ±  A)  -*  i  we  find 


»  rir  (X  t  j)  -r  0(A") 


Combining  these  results  gives 


»  ±  rirx  g(f^)  +  i  [r-rxA  g(fo)]'  +  O(A^) 


Thus 


2  C^Cf^)  =  [r-rxA  gCf^)]'  +  O(A^) 


Furthermore 


P  ±  f  [  1  -  x(l  -  2a  f^)  ]  4-  0(A) 


±  A  ‘  2 


which  gives 

♦  (X) 


sx 


1  +  i£  +  0(A) 


Substituting  these  results  into  (13)  converts  it  to  the 


perfect  differential  equation: 


s  X  f  ' 
o 


U  ^  F  J 


Integrating  we  obtain  (relevant  to  A  «  1) 

1 


f^(x) 


s  9 

X  S  *  ^ 

^  ^  ® 


(23b) 


(230 


(23d) 


(23e) 


(23f) 


(23g) 


following 


(24) 


(25) 


14 


ii 


where 


(25a) 


These  limiting  solutions  will  come  into  play  in  constructing  a  solution 
to  (13). 


15 


III. 


A.  Integration  Technique 

To  convert  (13)  to  a  more  tractable  form  we  write  (deleting  the 
subscript  on  f^) 

2^  C^(f)  -  g(F)]'  +  K(x)  (26) 

Note  that  the  lead  terms  in  this  approximation  maintains  the  property 
(21a)  and  is  consistent  with  the  finding  given  by  (23e) .  The  trial 
function  ')<(x)  which  also  contains  arbitrary  parameters  is  constructed 
in  a  manner  which  effectively  minimizes  the  error  term  k.  Having  thus 
constructed  a  solution  to  (13)  is  then  obtained  in  which  k  is 
neglected. 

B.  Approximate  Solution 

We  now  consider  the  solution  to  the  complete  equation  (13)  with 
the  approximation  (26)  and  k  =>  o. 

(♦f)'  +  [X^  g(f)]'  +  ['<'  g(f)]'  -  0  (27) 

where,  we  recall,  is  an  arbitrary  function  to  be  determined  and  ♦  as 
given  by  (14)  depends  on  f.  Integrating  (27)  gives 

(♦  +  x^  +  v)^)f'  +  (x^  +  '(»)f(l  -  af)  =  (28) 

where  is  a  constant.  This  equation  may  be  further  integrated 

provided  ♦  is  weakly  dependent  on  f.  This  is  the  case  when  a  f  «  n^^. 
In  the  classical  domain  where  a  »  0,  this  condition  is  trivially 
satisfied.  In  the  quantum  domain  (a  *  1)  this  condition  become  f  << 
n^^  which  is  compatible  with  the  property  f  <  1  and  the  fact  that  the 
distribution  n^^  relevant  to  bosons,  may  be  greater  than  1. 
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0,  we  obtain 


Assuming  this  constraint  is  satisfied  and  selting  « 

f  (X)  -  - - - - — - -  (29) 

X  ♦  (y) 

a  -f  B  exp  [x  -  /  [ - 2 - - - ]  dy 

o  +  y  "*■  '^(y) 

where  ♦q  equals  ^  with  a  =  o. 

C.  Full  Solution 

With  the  f^(x)  as  given  by  (29)  at  hand  we  now  seek  to  determine 

the  solution  to  (13),  f(x).  To  these  ends  the  LHS  of  (13)  is  rewritten 

L(f)  H  (♦f)'  +  [x^  g(f)]*  +  I  c.(f)  (30) 

±  * 

and  we  consider  the  functional 

00 

r[f]  s  /  [L(f) dx  ^  0  (31) 

o 

where  the  equality  corresponds  to  the  solution  to  (13).  The  final  form 
of  f(x)  is  determined  by  minimizing  r[f]  or,  equivalently,  by  finding  f 
such  that 

5  rcf]  =  0  (32) 

where  5  is  a  variation  with  respect  to  f. 

Motivated  by  the  results  (23  c,  e)  relevant  to  case  c,  (A  -*  0)  ,  we 

set 

=  b  +  (A  -  b)  f_  (33) 

which  at  b  «  A/2,  and  recalling  (26),  returns  the  asymptotic  form 
(23e)  ,  i.e. ,  ^  »  rirAx. 

Numerical  attempts  at  solving  (32)  at  A  <<  1  strongly  suggest  that 
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the  solution  be  separated  into  two  components  corresponding  to  x  <  A 


and  X  >  A . 
f(x)  - 


The  resulting  form  is  given  by 
c  f^(x)  ^  (1  -  c)  [f^(A)  +  f^'(A)(x  -  A)],  X  <  A' 

.f^(x)  ,  X  >  A 


(34) 


which  is  seen  to  be  continuous  at  x  «  A .  With  (33)  and  (34)  we  find 
that  the  functional,  r[f],  becomes  a  function  of  the  parameters  (b,c) . 
The  conditions  (32)  then  becomes 

ar  ar 

ab  *  dc  *  ^  (^5) 

Solutions  to  these  equations  determine  b  and  c  as  functions  of  the 
parameters  (T,  A,  s) . 

As  we  wish  to  compare  final  results  with  experiment,  this 
procedure  was  carried  out  for  n-type  Silicon  at  300°K  with  A  <<  1. 
Resulting  values  of  b  and  c  vs  electric  field,  s,  are  shown  in  Fig.  3. 

Approximate  analytic  expessions  describing  these  values  are  given 
by 

b(«)  =  0(s-48O)  logj^g  9  (480  -  s)  (36a) 


c(s)  =  I  (1  -  tanh  [logj^Q  (36b)  J 

Characteristic  behavior  of  these  curves  at  small  and  large  s  is  due  to 

the  following.  From  (14)  we,  note  that  ♦(x)  o>C  s.  Thus  at  small  s,  J 

recalling  (29)  we  see  that  f^(x)  grows  insensitive  to  4<(x),  and 
therefore  to  b.  This  property  was  confirmed  in  numerical  analysis 
where  it  was  observed  that  a^f/db^  2:  0  in  the  domain  s  i  100.  ■ 

Concerning  c(s) ,  it  was  shown  that  the  solution  to  the  principal 
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equation  (13)  goes  to  the  Fermi-Dirac  distribution  as  s  -♦  0.  Thus  with 

reference  to  (34)  we  see  that  this  behavior  corresponds  to  c  -»  i.  For 

large  s,  one  expects  the  component  f(x  >  A)  to  dominate  in  (34).  As 

2  2 

this  component  does  not  contain  c,  again  one  finds  d  V/dc  =:  0  in  this 
domain  (i.e.,  s  ^  lOO) . 

With  these  values  of  b(s)  and  c(s)  at  hand  the  distribution  f(x) 
as  given  by  (34)  is  plotted  for  various  values  of  s  in  Fig.  4.  From 
these  curves  it  is  evident  that  f(x  >  40)  (corresponding  to  A  >  leV)  is 
vanishingly  small.  This  is  consistent  with  our  omission  of 
non-equivalent  intervalley  scattering. 

0.  Drift  Velocity 

With  the  normalization  of  f(k)  as  given  by  (3)  we  write 

c  ”  (2ir) 


Converting  to  integration  over  x  and  4,  and  recalling  (16)  we  obtain 


Xd 


4  uVs  E.  I'"  x^  ^ 

'  V37  ^o  X  +  +  f_  dx 


dx 


(38) 


In  obtaining  (38)  we  employed  the  tensor  relation 

/  ick  dk  -  i  /  dk 


(38a) 


(recall  4  »  E;k;. 

Integrating  (38)  by  parts  gives  the  final  form  (with  S  =  +  f_) 


where,  we  recall,  f(x)  =  fQ(x)  is  given  by  (34).  A  plot  of  (39)  for 
n-type  Silicon  at  300°K  and  A  <<  l,  as  a  function  of  electric  field  E, 
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is  given  in  Fig.  5.  A  compilation  of  experimental  results^^  is  drawn 
on  the  same  graph,  it  should  be  emphasized  that  this  agreement  with 
experiment  in  no  way  depends  on  the  fitting  of  free  parameters. 
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As  a  self  consistent  check  of  the  solution  (34)  we  re-evaluated 
the  solution  numerically  employing  a  Monte  Carlo  technique  such  as 
described  in  reference  4.  Numerical  values  of  parameters  used  in  this 
analysis  are  given  in  Table  l.  The  analysis  addresses  the  steady-state 
solution  of  (1)  with  a  »  0  and  scattering  rates  described  in  Appendix 
A.  In  this  procedure  we  obtain  f(e,u).  The  t  *  0  component  of  (5a)  is 


then  given  by 


A  plot  of  f^(i)  so  obtained  for  the  specific  case  s  •  100  is  compared 
with  our  analytic  result  (34)  and  is  shown  in  Fig.  6.  The  discreteness 
of  the  Monte  Carlo  result  is  due  to  the  finite  mesh  in  energy  whereas 


the  fluctuations  of  the  results  are  due  to  the  statistical  nature  of 


the  numerical  procedure.  The  Monte  Carlo  curve  shown  corresponds  to  a 
sample  of  10,000  collisions. 
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V. 


A  nonlinear  difference-differential  equation  for  the  distribution 
of  charge  carriers  in  a  semiconducter  including  strain  acoustic  and 
equivalent  intervalley  scattering  was  derived.  Introducing  an 
approximate  differential  form  for  one  term  in  this  equation  permitted 
exact  solution.  This  solution  was  then  used  in  tandem  with  a  numerical 
analysis  to  obtain  an  approximate  solution  to  the  original  equation  in 
the  near-classical  domain.  The  distribution  so  obtained  was  applied  in 
evaluation  of  drift  velocity  for  n-type  Silicon  at  room  temperature. 
Results  were  found  to  be  in  very  good  agreement  with  measured  values 

5 

for  electric  fields  up  to  10  V/cm. 
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^pp^nqix  ft 

Scattering  Rates 

In  this  appendix  we  calculate  the  scattering  rates  that 

appear  in  our  starting  equation  (2)  and  the  net  rates  given 

by  (4)  . 

For  SA  interactions  we  write 

=kk“’  =  '"sa"3'  »  J  »  5I  «(«■  -  J  uq)  (All 


Where  upper  and  lower  signs  refer,  respectively,  to  phonon  absorption 
and  emission.  Furthermore 
±  q  -  k*  -  k 


represents  phonon  wavevector  and 

1 


ns;^(q)  = 


exp(^)  -  1 


-ftuq 


1 

2 


(A2) 


where  we  have  recalled  the  approximation  (8) .  For  IV  interactions  we 
write 


(IV±)  ^  2»  .D  rt 


'kk' 


(A3) 


For  the  net  scattering  rates  (4),  with  (3b)  we  write 


(A4) 


Inserting  the  preceeding  expressions  and  integrating  we  obtain 

x(SA±) 


k 


"^1  1  V  (2k±2mu/fi)^  -  1  (2kf2mu/r;)^. 

>  l-M.t  2  2  3 


.^2^„  k  '•fiu 
4irn  pu 


(A5) 
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and 


I 


x(IVi) 


-  V2m  ^TT- 

27rtt^puj^  -ft 


“IV®(*  ±  -Rujy)  (n^^  ^  i  ^  J) 


where  is  given  by  (9a) 


In  this  appendix  we  wish  to  obtain  the  intervalley  collision 

A 

integral  given  by  (9) .  Combining  (2)  and  (A3)  together  with  the 
limiting  process  (3a)  gives 


[f'(l  -  af)  (njy  -"1*7)  -  f(l  -  af)  (n^.^  +  i  +  i)  ] 


1-1. 


Now  we  note 


—  vn  dl'dM’d# 
2fi^ 


where  with  (5b) 


u '  *  k* • E 


and  ♦  is  an  azimuthal  angle  about  the  E  axis.  Furthermore 


f  s  f^d')  + 


with 


/dM*  =  2  ,  fdu'u'  "  0 

and  recalling  given  by  (9a),  integration  over  returns  (9) 
in  particular  that  the  difference  terms  in  (9)  stem  from  the 
function  in  (Bl) . 


Note 

delta 
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1. 


( 


aura  Captions 

Diagram  Illustrating  electron>phonon  interactions  relevant  to 
Silicon.  Notation  is  defined  in  the  text. 

2.  Scattering  rates  in  n-type  Silicon  at  300°K.  (a)  Strain-acoustic 

(SA)  and  g  intervalley  scattering.  (b)  f  intervalley  scattering. 
Absorption  and  emission  are  referred  to  as  (+/**)  respectively 
whereas  longitudinal  and  transverse  are  referred  to  as  L  and  T 
respectively.  Each  curve  for  corresponds  to  a  distinct  value 

of  a  whereas  k  is  the  wavevector  corresponding  to  energy  I . 

3.  (a)  Parameter  b  vs  log^^s,  as  it  appears  in  (33). 

(b)  Parameter  c  vs  log^^^s,  as  it  appears  in  (34). 

4.  Plots  of  the  distribution  f^(x)/A  for  various  values  of  electric 

field  s  in  the  limit  A  «  1.  The  optical  frequency  A  »  2.4  is 

shown . 

5.  Derived  values  of  drift  velocity  (S-L) ,  compared  to  experimental 
values  (Jacoboni  et  al^^) . 

6.  Monte  Carlo  (•)  and  analytic  ( — )  distribution  function  vs 
charge-carrier  energy  for  n-type  Silicon  at  300°K,  s  =  100  and  A 
<<  1. 
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avt. 


SA+ 


♦  t  \  * 

O  <<t<< 
-J 


39S  ^  <  jY )  0‘60| 


LA- 


62 


0  0 


( s/ujo !  Pa)  0>6o| 


V/cm) 


0.12 


(Aa);? 


